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Abstract 



We present a unified treatment of tliree cases of quasi-exactly solvable 
problems, namely, charged particle moving in Coulomb and magnetic fields, 
for both the Schrodinger and the Klein-Gordon case, and the relative motion 
of two charged particles in an external oscillator potential. We show that all 
these cases are reducible to the same basic equation, which is quasi-exactly 
solvable owing to the existence of a hidden s/2 algebraic structure. A system- 
atic and unified algebraic solution to the basic equation using the method of 
factorization is given. Analytic expressions of the energies and the allowed 
frequencies for the three cases are given in terms of the roots of one and the 
same set of Bethe ansatz equations. 



PACS: 03.65. Ge, 03.65.Fd 



1 



I. INTRODUCTION 



It is well known that exact solutions are hard to come by in physics (in fact, in all sci- 
ences). Many exactly solvable examples presented in textbooks of physics are only rare cases. 
More often than not they serve only as paradigms to illustrate the fundamental principles 
in their respective fields. For real problems approximation methods are indispensable. 

Recently, it was found that for certain quantum-mechanical problems analytical solutions 
are possible, but only for parts of the energy spectra and for particular values of the funda- 
mental parameters of the problems. First it was realized that the problem of two electrons 
moving in an external oscillator potemtial was of such a class Later, it was discovered 
that two-dimensional Schrodinger equations of electron moving in an attractive/repulsive 
Coulomb field and a homogeneous magnetic field also share the same characteristics [ 0- 
P]]. More recently, the latter problems were extended to the two-dimensional Klein-Gordon 
case and the Dirac case [0. 

The essential features shared by all these examples are as follows. The differential equa- 
tions (Schrodinger, Klein-Gordon, and Dirac) are solved according to the standard proce- 
dure. One first separates out the asymptotic behaviors of the system. One then obtains 
an equation for the part which can be expanded as a power series of the basic variable. It 
is at this point that deviation from the standard exactly solvable cases appears: instead 
of the two-step recursion relations for the coefficients of the power series so often encoun- 
tered in exactly solvable problems, one gets three-step recursion relations. The complexity 
of the recursion relations does not allow one to do anything to guarantee normalizability 
of the eigenfunctions. However, one can impose a sufficient condition for normalizability 
by terminating the series at a certain order of power of the variable; i.e. by choosing a 
polynomial. By so doing one could obtain exact solutions to the original problem, but only 
for certain energies and for special values of the parameters of the problem. For the works 
mentioned above, these parameters are the frequency of the oscillator potential and the 
external magnetic fields. 
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Soon after it was realized that the above quantum-mechanical problems are just 
examples of the so-called quasi-exactly solvable models, recently discovered by physicists 



and mathematicians [ [Q- [T^]. This is a class of quantum-mechanical problems for which 
several eigenstates can be found explicitly. The reason for such quasi-exactly solvability is 
usually the existence of a hidden Lie-algebraic structure [ |TD|- 10]. More precisely, quasi- 



exactly solvable Hamiltonian can be reduced to a quadratic combination of the generators 
of a Lie group with finite-dimensional representations. 

In this paper we would like to show that three of the four problems mentioned in the 
second paragraph, namely. A) charged particle moving in Coulomb and magnetic fields 
(Schrodinger case); B) charged particle in Coulomb and magnetic fields (Klein-Gordon case); 
and C) relative motion of two charged particles in an external oscillator potential, can be 
given a unified treatment. We shall show that all these cases are simply variations of the 
same basic equation [eq.(0) below], which is quasi-exactly solvable owing to the existence 
of a hidden s/2 algebraic structure. This algebraic structure was first realized by Turbiner 
for the case of two electrons in an oscillator potential We shall give a systematic and 
unified algebraic solution to the basic equation using the method of factorization presented 
in [[7| for the case A. Our method allows one to find the analytic expressions of the energies 
and the allowed frequencies once and for all in terms of the roots of a set of Bethe ansatz 
equations. This is in sharp contrast to the method of solving recursion relations, which 
must be performed for each and every order of the polynomial part in order to get these 
expressions. Our treatment also reveals that the eigenenergies and the allowed frequencies 
in all these cases are given by the roots of the same set of Bethe ansatz equations. This also 
makes explicit the close connection between the three cases. 

We will define the three problems in Sect. II. The basic equation underlying these three 
cases is then solved in Sect. Ill by the method of factorization. In Sect. IV we show how 
previous results, obtained for the three cases by solving recursion relations, can be easily 
reproduced from the solutions obtained in Sect. III. The Lie-algebraic structure underlying 
the basic equation is discussed in Sect.V. Sect. VI then concludes the paper. 



II. THE THREE CASES 



In this section we shall give a brief description of the three cases of charged particles 
moving in external fields which we will consider in the rest of the paper. Following previous 
works, we adopt the atomic units h = m = e = lm the CGS system. 

A. Electron in Coulomb and magnetic fields: Schrodinger case 

This general case was considered in [ [Q- ||5[, [0]]. The Hamiltonian of a planar electron 
in a Coulomb field and a constant magnetic field B = Bz {B > 0) along the z direction is 

where c is the speed of light, Z (positive or negative) is the charge of the source of the 
Coulomb field, and the vector potential A is A = |B x r in the symmetric gauge. 
An ansatz of the eigenfunction in the polar cordinates (r, 6) is 

^{r,t) = ^exp{ime -iEt) , m = 0, ±1, ±2, . . . (2) 

Here m is the angular quantum number, and E the energy. The radial wave function u{r) 
satisfies the radial Schrodinger equation 

'1 d"^ 1 / o 1\ 1 loo Z 



u{r) = (3) 



where u;^ = i?/2c is the Larmor frequency . 



B. Electron in Coulomb and magnetic fields: Klein-Gordon case 

In the above problem is extended to the Klein- Gordon case, assuming the same 
ansatz of the wave function as in (0). Now the radial wave function u{r) obeys the following 
equation: 



1 1 f , Z^ 1\ 1 1 ^ , EZ c"^ 



u{r) = . (4) 
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But now, as noted in [H, the quantum number m must satisfy the relation 



m ^ > 



(5) 



in order for the solutions to make sense. This relation forbids the existence of the s-states 
(m = 0). 

C. Relative motion of two electrons in an external oscillator potential 



In [0 the author considered the problem of three-dimensional Schrodinger equation for 
two electrons (interacting with Coulomb potential) in an external harmonic-oscillator po- 
tential with frequency uJext- The Hamiltonian is 

1 



11 11 

H = "2 Vi + 2^extri - 2^2 + 2^extT^2 



(6) 

|ri - ral 

The total wave function is factorizable into three parts which depend respectively only on the 
center of mass, the relative coordinates, and the spins of the electrons. The wave function 
of the center of mass coordinates satisfies the Schrodinger equation of a three-dimensional 
oscillator the solution of which is well known. The spin part dictates the parity of the wave 
function of the relative motion. The Schrodinger equation for the relative motion is 



^Y72 I ^, ,2„2 , 1 

"2 V. + 2'^rr + ^ 



e(f){r) 



(7) 



where r = ri — T2, oJr = c<Jext/2, and e' is one half of the eigenenergy of the relative motion (in 
the notation of 0). By assuming an ansatz of the wave function in the spherical coordinates 
of the form 



0(r) = ^l^^(f) , 
r 

where Yim are the spherical harmonics, we get from (|^ the following equation 
'Id? /(/ + 1) 1 1^0 1 



11 1 2 2 
--2 - t'rr 



+ e 



'u(r) = . 



(8) 



(9) 



2c/r2 2 r2 2"'" 2r 

We note here that if we change the sign of the 1/r term in the last equation, we get an 
equation that describes the relative motion in the oscillator potential of an electron and a 
positron. This case is included in our discussions. 



III. THE BASIC EQUATION AND THE METHOD OF FACTORIZATION 



It can be seen that equations (|), (| 
have the same basic form, namely: 

7(7-1) 1 



and (y), after an appropriate change of parameters, 



2c/r2 



1 2 2^/3^ 

-to r H ha 

2 r 



u{r) = . 



(10) 



Here /3,7 and u (7,0; > 0) are real parameters, and a is the eigenvalue of eg. (PDf) . Exphcit 
expressions of these parameters for the three cases mentioned in the Introduction will be 
given in the next section. That this equation is quasi-exactly solvable means that, given a 
fixed value of the parameter 7 and P (or u), the equation can only be solved exactly only 
for particular set of parameter u (or /?) and eigenvalue a. 

Now we make the following change of variables: x = v^2cJr and b = J2/uj[3. Then 



eq.(lO) becomes: 



d'' 7(7 - 1) 



X b a 
— + - + - 

A X UJ 



U[X] 



. 



(11) 



The values of a and b in eq.(|TT]) may be found by means of a method closely resembling 
the method of factorization in nonrelativistic quantum mechanics 0. We shall discuss this 
method briefly below. Let us assume 



u{x) = x'*' exp(— x^/4)Q(x) 



(12) 



where Q is a polynomial. As mentioned in the Introduction, the assumption that Q be a 
polynomial is only a sufficient, not necessary, condition for normalizability of the eigenfunc- 
tion u{x). Substituting (0) into (pUD , we have 



dx^ \ X ) dx \ X , 



Q{x) = 



(13) 



where e = a/cij — (7 + I/2). 

It is seen that the problem of finding spectrum for 
eigenvalues of the operator 



is equivalent to determining the 
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H 



27 \ d h 
--Vd^-'x 



(14) 



We want to factorize the operator in the form 



H = a^a + p, 



(15) 



where the quantum numbers p are related to the eigenvalues of (|T3D by p = e. The eigen- 
functions of the operator H at p = must satisfy the equation 



aip = 



(16) 



Suppose polynomial solutions exist for (ITSD, say Q = Yl {x — Xk), where Xk are the zeros 

k=l 

of Q, and n is the degree of Q (we mention here that the order n in this paper is equal to 
{n — 1) in [ [0]- 0] where x"'~^ is the highest order term in Q). Then the operator a must 
have the form 



d " 1 



k=l 



X Xj^ 



(17) 



and the operator a+ has the form 



^ 27 ^ 1 

+ x-2^ 

ox X f~, X — Xk 



(18) 



k=l 



Substituting ( p^Tf ) and (|T8D into (p!5|) and then comparing the result with (|T^), we obtain 



the following set of equations for the zeros Xk (the so-called Bethe ansatz equations |ll3|): 



as well as the two relations: 



Xfc - 2 2^ 



j^k -^i •'^k 



, k = 1, . . . ,n , 



(19) 



(20) 



k=l 



Summing all the n equations in (|T9D enables us to rewrite the first relation in (|^) as 



b=^Xk . 

k=l 



(21) 
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From p = e and the second equation in (^) , one gets 

e = n = a/^ - (^7 + . (22) 

For n = 1,2 the zeros Xk and the values of the parameter b for which solutions in terms 
of polynomial of the corresponding degrees exist can easily be found from (p!9| ) and ( |2TD in 
the form 



n = 1 , xi = ±a/27 , b = ±a/27 



n = 2 , Xi = —X2 = a/27 + 1 , 6 = 



x, = 2j/x2, a;2 = ±(l + V47 + l)/V2, 6 = ±^2(47 + 1) . (23) 

These expressions will be employed in the next section. 

For general values of n solving for the x^s from eq.(0) is difficult, and one must resort 
to numerical methods. But some properties of the solutions are known. First, from eq.(|l9D 
we see that, if {xk} is a set of solution to (p^D, then so is {—Xk}- This means, by eq. , 
that for every possible value of b, there is a corresponding negative value —b. Second, as we 
shall see later in Sect.V, the number of values of b for a fixed order n is n + 1. 

IV. SOLUTIONS TO THE THREE CASES 

We now apply the results of Sect. Ill to reproduce the results of previous works. The 
essential step is to solve the Bethe ansatz equations ([19|) for the roots x^^s for each order 
n. Then from eqs.(^) and (|22D we obtain the values of the allowed pair of frequency and 
energy. It will be apparent that it is the values of 6's that determine the energies and the 
frequencies u = 2/?^/ 6^. 

A. Electron in Coulomb and magnetic fields: Schrodinger case 

In this case, 7 = \m\ + 1/2, u = ujl, (3 = Z = ±|Z|, and a = E — muji- The upper 
(lower) sign in j3 corresponds to the case of attractive (repulsive) Coulomb interaction. This 
will be assumed throughout the rest of this subsection. We have 

8 



= — ^ , E = LJL{n + m+\m\ + l) . (24) 

These are the general expressions for the frequency (and hence the magnetic field) and the 
energy in terms of the values of b. 

The case of = 1 and n = 2 can be obtained easily. From (^Sf) and the definition of b 
one has 

n = l, o;z. = 2 ^ = o\ i ^ 3 + ^ + i 

l\m\ + 1 l\m\ + 1 

n = 2, u;l= . i^2= ■ ^ 4 + m+m ■ 25 

4|m| + 3 4|m| + 3 

The corresponding polynomials are 

Qi = X — Xi = X =F \/27 , 

2 



^2= n(^-^fc) =^'Ta;v/2(47+l) + 27 . (26) 
fc=i 

The wave functions are described by (|I^. For the repulsive Coulomb field the wave functions 
(for = 1,2) do not have nodes, i.e. the states described by them are ground states, while 
for the attractive Coulomb field the wave function for n = 1 has one node (first excited 
state) and the wave function for n = 2 has two nodes (second excited state). 

Let us mention here that we may consider a dual situation of the original problem: we 
may consider the magnetic field B (and thus cj^) as a fixed quantity, and eq.(p^ then give 
the allowed values of the energy and the charge Z. 



B. Electron in Coulomb and magnetic fields: Klein-Gordon case 

For definiteness, we consider positive energy solutions for the attractive Coulomb po- 
tential (Z > 0). This is the case considered in Negative energy solutions, and the 
case for repulsive Coulomb field can be treated in exactly the same way. In this case, 
7 = ^m^ - Z'^jc^ + 1/2, = cji, /5 = ZEjc^, and a = E^/2c'^ - c^/2 - muji- In order for 
the wave function to make sense, 7 has to be real. This implies that — Z'^ j(? > 0, which 
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forbids the existence of the m = states (the s states) in the Klein-Gordon noted 



m 



Using ul = 2/3 /b we get the allowed magnetic field as 



B = 2cul 



(27) 



and the corresponding energy E is obtained from eq. 



n 



a 



(7 + 1/2) 

62c6 1 



4^2 4^2 E^ 



1 + m+ J m? — Z^/c" 



(28) 



This equation leads to 



E' 



6V 



n + l + m + dm? — Z"^ /c 



-1 



(29) 



These are the most general expressions for the energy and the frequency. 

To compare with the results of = 1, 2 cases given in H, we substitue these values of n 



into dH) and 



For n = 1, we get 



E = c' 



B 



AZ^ [2 + m+ Jrn? - Z'^/c'^ 



c2 {2^m^ - Z^/c^ + 1 
AE^Z^ 



-1/2 



c3 (2^m2-ZVc2 + l) 



(30) 



and for n = 2, 



E = ^ 



B 



2Z2 (3 + m + ^w? - Z^/c 

c2 (4^m2- Z2/c2 + 3^ 
2E^Z^ 



-1/2 



(31) 



c3 (^4^m2-Z2/c2 + 3j 

These expressions are exactly the ones obtained in by solving recursion relations. 

For negative energy solutions, the energy is given by the negative roots of eq. (pOj) . The 
only difference is that the roots of the Bethe ansatz equations have opposite signs, in view of 
eq. (|2TD . This only changes the nodal structure of the wave functions. From the expression 
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P = ZE/c^, we note the equivalence between the positive (negative) energy solutions in 
the attractive Coulomb case and the negative (positive) energy solutions in the repulsive 
Coulomb case. 

Similar to case A, one may consider the dual situation in which the magnetic field is 
assumed fixed, and the Bethe ansatz equations instead give the values of the allowed pair 
of energy and the Coulomb charge. 



C. Relative motion of two electrons in an external oscillator potential 

In this case, •j = I + 1/2, uj = Ur, P = —1/2, and a = e' . We have the following general 
solutions: 

1 / 3\ 

^r = ^ , e' = ujr[n + l + - ) . (32) 



For the two simplest cases n = 1 and n = 2, we have 



1 , 2/ + 5 

n = 1 , ujr = -TT, — — , e 



4(/ + l)' 8(/ + l) 

-9 1 , 2/ + 7 

" ' " 4(4/ + 5) ' ' " 8(4/ + 5) 



(33) 



These are exactly the expressions given in p|. They are also the solutions (for n = 1, 2) for 
the case of an electron and a positron in the oscillator potential {(3 = +1/2). The nodal 
structures of the wave functions are the same as those described for case A. 



V. HIDDEN LIE- ALGEBRAIC STRUCTURE OF THE BASIC EQUATION 

The basic equation ([TD|), or its equivalent form (|T^, possesses an underlying Lie-algebraic 
structure that is responsible for its quasi-exactly solvability. 

In fact, Turbiner has indentified a s/2 structure for the case of two charged particles in 
an oscillator potential p. In view of the fact that all the previous cases considered in this 
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papers are related to the same basic equation (|IOD, one expects the same hidden structure 
to be present in all these cases. This is indeed the case, and it is sufficient to show that s/2 
algebra is in fact the underlying structure possessed by or (|13|). In this section we shall 
carry out Turbiner's analysis to eq.(0), with only slight modifications in the parameters to 
suit the general situation. Only the main ideas are given here, and we refer the reader to 
§ for details. 

Let us construct three generators in the following manner: 

T+ 2 ^ 

Jn = r — - nr , 
dr 

o_d_n 

" ~ dr 2 ' 

d , s 

^» = * ■ 

These generators realize the sli algebra: 

for any value of the parameter n. If n is a non-negative integer, then there exists for the s/2 
algebra a (n + l)-dimensional irreducible representation P„+i(r) = (l,r, r^, . . . ,r"). From 
this it is clear that any differential operator formed by taking polynomial of the generators 
(0) will have the space Vn+i as the finite-dimensional invariant subspace. This is the main 
idea underlying the quasi-exact ly- solvable operators [ [§- |]14| ]. 

Now consider the quasi-exactly-solvable operator which is quadratic in the J„'s: 

T2 = - J° J- + 2^ J+ - Q + 27) 4^ - un . (36) 



This operator belongs to the class VIII according to the classification given in |T0[ . In terms 
of r, T2 becomes 

T2 = -r— + 2 Lr^ -l)-r- 2w«r . (37) 
dr^ ^ ^ dr 

Let us now consider the eigenvalue problem 

T^Qir) = 2f3{n)Q{r) . (38) 
12 



This eigenvalue problem possesses n + 1 eigenvalues /3(n), and the corrsponding eigenfunc- 
tions are in the form of polynomial of the n-th power, while other eigenfunctions are non- 
polynomial which in general cannot be found in closed analytic form [^]. Let us first substi- 
tute the form ( pTf ) into eg. (pB]) , then divide the resulting equation by r/2, we arrive at 

' ( ^ ' ■ ^ ■ Q(r) = . (39) 



2 dr"^ \ r J dr ~^ ~^ r 

Finally, we divide (|39|) by lu and change the variable r to x = ^/2ujr. This leads to the 
following equation 

d^ f2j \ d I /2 , 1 , , 

' ' ' , 1^ , / 1 Q(2;) = . (40) 



dx"^ \ X J dx \ \ uo X 



This is exactly eq.(|T3D, provided that e = n and h = ^2/u!P{n). This means that eq.(|T3D is 
quasi-exactly solvable if e = n, which is exactly our relation in eq. (p^) , and that there are 



only n + 1 allowed values of 6 = ^J2/uP in eq.(|l^) [cf. eq. (p3D]. 

Translating back to the original three cases considered in this papepr, these results 
imply the following. In case C, /3 is a fixed parameter {f3 = ±1/2), hence the finite number 
{= n + 1) of the values of b implies the same number of the allowed frequency Uext of the 
external oscillator potential and the corresponding energy. This is the case found in [0] and 
presented here again from a new light. For case A {f3 = Z) and case B (/3 = ZE/c^), the 
above results mean that, at a fixed order n, there are exactly n + 1 allowed values of the 
pair of energy and magnetic field for a fixed Coulomb charge, or of the pair of energy and 
Coulomb charge for a fixed magnetic field. 

Furthermore, it has been shown ^ that there exist [{n + l)/2] positive eigenvalues and 
the same number of negative eigenvalues of the b (here [a] represents the integral part 
of a). In the general situation considered in this paper, positive (negative) values of b 
correspond to the attractive (repulsive) Coulomb field for positive energy solutions. For 
negative energy solutions, the sign of b is reversed for the two kinds of Coulomb field. 
Hence, our unified treatment together with the Lie-algebraic analysis of these cases give a 
very simple explanation as to why the number of the positive energy levels for a fixed order 
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n considered in [ 0- 0] are all equal to [(n + l)/2]. 

VI. CONCLUSIONS 

In this paper we have presented a unified treatment of three cases of quasi-exactly solvable 
problems, namely, charged particle moving in Coulomb and magnetic fields, for both the 
Schrodinger and the Klein-Gordon case, and the relative motion of two charged particles 
in an external oscillator potential. We show that all these cases are reducible to the same 
basic equation [eq. (p!0D ], which is quasi-exactly solvable owing to the existence of a hidden 
SI2 algebraic structure. A systematic and unified algebraic solution to the basic equation 
using the method of factorization is given. Our method allows one to express the analytic 
expressions of the energies and the allowed frequencies once and for all in terms of the roots 
of a set of Bethe ansatz equations. Our treatment also reveals that the eigenenergies and 
the allowed frequencies in these cases are all given by the roots of the same set of Bethe 
ansatz equations. 
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